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1 Abstract 

In this paper, a general equation for the evolution of an axisymmetric magnetic field in a Hall 
plasma is derived, with an integral similar to the Grad-Shafranov equation. Special solutions 
arising from curvature - whistler drift modes that propagate along the electron drift as a Burger's 
shock, and nonlinear periodic and soliton-like solutions to the generalized Grad-Shafranov integral 



- are analyzed. We derive analytical and numerical solutions in a classical electron-ion Hall plasma, 
in which electrons and ions are the only species in the plasmas. Results may then be applied to 
the following low-ionized astrophysical plasmas: in protostellar disks, in which the ions may be 



coupled to the motion of gases; and in molecular clouds and protostellar jets, in which the much 
^ . heavier charged dust in a dusty Hall plasma may be collisionally coupled to the gas. 
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I. INTRODUCTION 



Hall physics is relevant in plasmas in which one or more species of a plasma are non- 
magnetized, such as in laboratory plasmas on time scales much shorter than the inverse ion 
cyclotron frequency jl-5]. Within the past thirty years there have been numerous theoreti- 
cal, experimental, and numerical studies of Hall effects in plasmas. Experimental, numerical, 
and theoretical research on Hall plasmas have focused on planar rather than more compli- 
cated geometries. 

Whist.er drift .odes propagate a8 a shoe, wave via density grad.ents in a Hall p.a 8m a Q. 
However, the same drift modes, which propagate as shocks, and whistler-like modes with 
periodicity, may also be seen when considering curvature effects in the absence of density 
gradients. Whistler-like modes with curvature effects have been shown to be important in 
magnetic reconnection in the earth, magnetosphere flQ, and Hal. physics is important 
in characterizing reconnection |9j. Linear studies of Hall physics in a low- ionized, high- 



density protostellar disk 10|, [ll| (i.e. cylindrical geometry) have shown that the whistler 



frequency is also an important parameter in the stability of the rotating disk, and compre- 



hensive numerical studies 
Rudakov and others 



14 



13) have borne this out. 



161 ] have analyzed whistler and whistler drift modes with planar 



geometry and inhomogeneities in astrophysical dusty plasmas. This paper continues work 
done by Rudakov l_7j] and analyzes the nonlinear axisymmetric structure of whistler and 
whistler drift modes. For simplicity, we construct Hall magnetohydrodynamic (MHD) solu- 
tions for an electron-ion plasma, but initially demonstrate their applicability to low-ionized 
astrophysical plasmas. 

The organization of this paper is as follows. In Sec. (JTTJ) , we review appropriate dimen- 
sional normalizations (length, time, speed, magnetic field) for an electron-ion Hall plasma. 
In subsequent subsections we describe changes in physics, and corresponding changes in 
scalings, that expand the applicability of Hall MHD to low-ionized electron-ion-gas plas- 
mas (protostellar disks) and to dusty plasmas (protostellar jets ). In Sec. ( 1IIII) we derive 
the generalized Grad-Shafranov equation, first introduced in [17|, as an interesting class of 
solutions to an axisymmetric Hall plasma. In Sec. (llVj) we discuss interesting new solutions 



to the GGS beyond those described in 



17j ]. The concluding section, Sec. (jV]), summarizes 



this paper's results and illuminates an alternative explanation of what has been though to 



2 



be magnetic bubbles 



18[ moving along protostellar jets. 



II. HALL PHYSICS AND LENGTH SCALES 

In a standard electron-ion Hall plasma, magnetic fields are tied to the electron, rather 
than fluid, velocity. As a result, unlike in MHD, we may get penetration of the magnetic 
field into the bulk medium. Furthermore, for this traditional Hall plasma, one gets these 
scalings 

• The frequency of Hall phenomena are normalized in units of u C i = ZicB^/ (rriic), the 
ion cyclotron frequency. Z^e is the charge of the ion, Bq is the magnetic field strength, 
rrii is the proton mass, and c is the speed of light. 

• Length scales are in units of di = c/u P i, the ion inertial skin depth, where = 
47re 2 rii/mi is the ion plasma frequency, rij is the ion number density. 

• Velocities of Hall magnetic structures are in units of va, where va = Bo/^Annrrii is 
the Alfven velocity. 

• Magnetic diffusion is given by the following expression, 

V m = (1) 

v ei is the electron-proton collision frequency and m e is the electron mass. 



These scalings are also summarized in Eq. (fT2|) . Sec. fill Ap and fill Bp describe the low- 



ionized astrophysical plasmas in which Hall MHD phenomena may be observed. Different 
physics results in different dimensional scalings of Hall MHD phenomena. Most notably, 



in these low-ionized plasmas, magnetic diffusion is achieved throug 
collisionality, rather than electron- ion collisionality (see, e.g., 111. Il4 



i ion-gas and dust-gas 



19| for a review). 



A. Low Ionized Hall Plasmas 

In low ionized plasmas, charged particles make up only a minority of the number and mass 
density of the total medium. In general the electric field consists of inductive, Hall, Ohmic, 
and ambipolar components. In order that the low-ionized plasma may be approximated by 



a resistive Hall MHD electric field (Ohmic and Hall electric fields) with negligible ambipolar 
electric fields, the ions must remain collisionally coupled to the gas while electrons must 
remain collisionally uncoupled. Using estimates of those physical regimes in a low-ionized 



plasma where either the inductive, Hall, ambipolar, or Ohmic electric fields dominate [11], 
resistive Hall MHD plasma dynamics may occur where v eg v ig 3> u C iU ce and u ce 3> v eg . v ig 
and v eg are the ion-gas and electron-gas collisional frequencies, and u ce is electron cyclotron 
frequency. 

In low-ionized plasmas whose electric fields are resistive Hall MHD in nature, the following 
substitutions to an electron-ion Hall plasma are made. 

• We replace v ei with v eg in Eq. §1$. 

• The ions are coupled to the gas, so that the Alfven velocity is calculated over the den- 



sity of gas particles va 9 — B / ^4:7im g n g , where m g and n g are the averaged molecular 
mass and number density of gas particles, respectively. 

Since Hall plasma time scales are scaled with respect to the ion cyclotron frequency, 
the Hall length scale (the largest length scale over which we can observe Hall effects) 
is given by this relation: 



VAg c / m g y- I n g 



u ci uj p i \rrii j \n 



1/2 / \ 1/2 

(2) 



The collision rates for electrons and ions are given in 20j, |2l|, where A is the molecu- 
lar/atomic mass of the principal ion: 



Vig =n g (av) i 
=2.6 x 10 
=n g (av) eg 



><j 

-9 1 11 \ „-l 



A Y 1/2 - ( 3 ) 



lamu 



=8.28 x KT 10 — s" 1 



IK J 



1/2 (4) 



In a typical protostellar disk, n g ~ 10 12 — 10 13 cm 3 22j, magnetic fields are of order 10 



- 100 mG (see 



for scaling of magnetic fields in molecular clouds), so that we have a 



Hall plasma. Most of the ions in low-ionized astrophysical plasmas are alkali ions 22|, |24| , 
giving an ion mass mj ~ 30. For T ~ 30K, gas density n g = 10 13 cm -3 , assuming a standard 



cosmic-ray ionization £ = 10 s [25|, ion density rii ~ 1 cm 



10 



cm 



26|. Here z/ ifl ~ 1.4 xlO 5 
For B = 100 mG, u C i = 30 s" 1 <C v™. Furthermore, since v e i = 2.6 x 10 4 s _1 and 



uj ce = 1.8 x 10 6 s -1 , t'eg^g 3> Wci^ce and we may completely neglect the ambipolar electric 
field in the induction equation. The upper threshold for Hall effects L ~ 5 AU (7.5 x 10 13 
cm), comparable to the scale of temperature and density gradients in the outer parts of a 
protostellar disk. 



B. Physical Regimes of Dusty Plasmas 

Approximately 1% of the interstellar medium in our galaxy is composed of dust . 



Rudakov [14| showed how to apply the results of electron-ion Hall plasmas to dusty plas- 
mas. It is important to note that although dust may not be a primary charge carrier in 
a low-ionized astrophysical plasma, their large masses, small charges, and relatively small 
number densities opens up a length and time scale over which the dust remains unmagne- 
tized (whether through long dust gyroperiod or through collisional coupling with interstellar 
gas) while the faster and smaller-scale electrons and ions remain fixed to the magnetic field 
- hence dusty Hall MHD physics. The following changes to the standard electron-ion Hall 
MHD scaling must be made, 

rii ->■ n d 

eZi ez d 

(5) 

Vi -> v d 



Aire 2 ni [{z d n d /n,j) + (u ig /u ci ) ] 

n d and m d are the number density and mass of dust particles (of a single mass and size), z d 
is charge per dust particle, and v d is the velocity of the dust. In the limit of low temperature 
olasmas (T < 100 K), z d — — 1. Dusty plasma behaves as a Hall plasma if n d /rii ^> v ig jbj ci 



14l | . The dust inertial scale, the lower limit over which we may apply Hall physics, is 
(M d c 2 / ATre 2 n d z^) 1 ^ 2 = c/u pd , where u pd is the dust plasma frequency. 

In the limit that the Lorentz force on a dust particle (ev d B/c) is smaller than the drag 
force acting on the dust (a 2 v g)t h m g n gV d) where a is the size of a dust grain): 

< a 2 v g n g (6) 



eB 2 



m g c 



Then the dust is collisionally coupled to the gas. In this limit, the dusty plasma Hall 
length scale c/u p d increases by a factor (mg/rrid) 1 ^ 2 (rig/rid) 1 ^ 2 , similar to the length scale 
of the low-ionized gas in Eq. (j2J). m g and n g are the masses and number densities of gas, 
respectively 

Dense molecular clouds have a hydrogen gas density of n# ~ 10 3 — 10 6 cm -3 (see, e.g. 



19l . 127] ). Given a standard cosmic ray ionization £ 
K, this implies an ion density rii ~ 10~ 3 



10 2 cm 3 



0- 17 s- 1 



23 







25| at temperature T = 30 



26j | . For these gas densities, the 



271 ] . For a representative molecular 



magnetic fields are on the order of 10~ 5 — 10~ 3 G 
cloud density rin = 10 3 cm -3 at T = 10 K with representative grains of size a = 0.03 /im, we 
get densities [26| n e 



10" 



cm 



rid 



10 cm , and magnetic fields B = 10 G. 
Assuming densities of dust at 3 gm/cm 3 , this gives an average dust mass rrid = 0.5 x 10~ 15 

1 /2 

gm. Furthermore, the length scale of Hall phenomena is (M^c 2 /A r Ke 2 ndZ d \) = c/oo p d ~ 1000 
AU. Although the density and temperature model inside protostellar disks is still open to 
question, we have somewhat better estimates of the density of gas at the edges of protostellar 
disks and in protostellar jets. We are taking the grain radius of a = 1 fim that scatters visible 
light efficiently Since we have marginal extinction due to the scattering of light by the dust, 
RndO 2 ~ 1 and so ~ 3 x 10" 7 cm" 3 . The length scale for Hall phenomena calculated over 
this mass and density is c/u p d ~ 1000 AU. However, this length scale is smaller than that 
calculated because the Hall physics is dominated by smaller dust particles. The usual dust 
size distribution function is n d (a) oc a -3,5 j^, and the length scale c/oo p d changes as a 3 . 



III. HALL DYNAMIC EQUATIONS 

We consider the electron-ion Hall plasma, using a cylindrical geometry and fluid equation 
for electrons. The Hall dynamic equations in this form were first written in jl^]. Here, we 
begin with the Hall frozen- in equations with magnetic diffusion expression T> m in Eq. ([I]), 
and letting ions remain motionless. 

n lL (—) + nv e ' V ( — J = (B • V) v e - VV m x ( V x B) + £> m V 2 B 

(7) 



dt \ n J \ n 



„ „ 47r 4nen 
V x B = —J = v, 

c c 



We consider an axisymmetric magnetic field in cylindrical geometry. Only two variables 
completely describe the magnetic field: 1) the toroidal magnetic field B = B^; and 2) 
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the toroidal vector potential A = A^. This implies an expression for the three-component 
(poloidal and toroidal magnetic field): 

B = Be^ + V x (Ae ) 

dA 

B r = [V x (Ae t )\ T = —g- 
B z = [V x (Av)], = i^tt 
-6</> = B 

Where B r and B z are the radial and vertical components of the poloidal magnetic field, and 
B^ is the toroidal component. 

The electron velocity in terms of the magnetic field is shown here, 

c 



-V x B 



Alien 

< = ~t~~ — t v x ( Be *)l = ir- d 4- 

Anen r Anen dz 

z _ c c djrB) 

e Alien z Airenr dr 

s c r „ .„ . , ... c (Id ( dA\ A d 2 A 

^ = - 7— - [V x (V x (Ae ))] , = — - -- r— - - + 



(9) 



An en ^ Alien \r dr \ dr J r 2 dz 1 

Where u£, v*, and vf are the radial, vertical, and azimuthal components of the electron 
flow velocity, respectively With Eq. (Q, the continuity equation for the electron fluid, and 
Maxwell's equations, the plasma number density is steady-state dn/dt = 0. 



A. Hall Dynamic Equations 



Let V m be a constant. Employing Eq. (jSJ) and Eq. (J2J) we recover these equations, 
dB _ fd(rB) d d(rB) d \ B fl d f dB\ B t d 2 B 



dt nr 2 \ dz dr dr dz J \r dr \ dr J r 2 dz 2 



(10) 



dt \ dz dr dr dz J nr \r dr \ dr J r 2 dz 

d{rA) d d{rA) d \ 1 fl d f dA\ A d 2 A 
dz dr dr dz J nr \r dr \ dr J r 2 dz 2 
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if we make these scalings of time, space, magnetic field, density, velocity, and magnetic 
diffusion D m as shown below, 

t ->■ tu ci (r, z) -)■ (r, z)u)pi/ c 

B ->• A ->• A/A = Auv/ (cB ) 

n^n/riQ ui ci = eB /miC 

= 47re 2 rao/^i f — > v/va 
D m -> z/ ei /W f a = -B / v/47rn mi 
cu ce = eB /m e c 

Bq and no are the maximum toroidal magnetic field (or magnetic field strength) and number 
density, respectively. u ce , uj ci) u P i, va, are the electron cyclotron frequency, ion cyclotron 
frequency, ion plasma frequency, and Alfven velocity at that point, respecively. 



B. Generalized Grad-Shafranov (GGS) Integral 

We show an important integral of the Hall dynamic equations, following much of the 
derivation shown in [l^]. We assume propagating solutions A = A(r,z — ut) and B = 
B(r,z — ut), where u is a constant normalized velocity, and a radial- dependent number 
density n = n(r). 

Consider force balance equation in some moving frame in which the electric E' and 

magnetic B' fields are stationary in time. Electron gas pressure and resistivity are ignored. 

The electron velocity in a dimensionless form, using Eq. (jSJ), is given by: 

f f d(rB) IdB 1/19/ dA\ A d 2 A\ , . 

v e = V x B = K —^e z + -— e r + - -— r— - - + — - e^l3 

n nr or n oz n \r or \ or J r A oz z J 

Begin with the dimensionless frozen-in law for the electron fluid neglecting electron pressure: 

E + v e x B = (14) 

The propagating solution is moving with sufficiently small velocity to ignore relativisitic 
effects so that: 

E' = E + uxB, B' = B (15) 

In this comoving frame the electric field E' = 0, since electrostatic potentials are neglected. 
In the laboratory frame the electric field is given by: 

E = -uxB (16) 



As a result, Eq. ()14p can be rewritten as: 

(v e - u) x B = (17) 

A general case reproduces the generalized Grad-Shafranov (GGS) integrals with force bal- 
ance in r, z, and <fi directions in Eq. ( I17p . Sufficiently, force balance in the radial direction 
implies: 

1 d ( dA\ A d 2 A\ ( ld(rB) 



r dr y dr ) r 2 dz 2 ) z ~^\~^rdr) ^ ^ 

Combined with the fact that the solution consists of magnetic surfaces at constant poloidal 

flux \I/ = rA, and that fact that dH/dr = dH/d^ x d^/dr, imply the GGS integrals. 

d 2 ^ dH[ u .^ [ r . , w ,,\ , inN 
r ¥r{r^) + ^ = 'd^{ H ^- U l n(r)r " r J (19) 

rB = H{^)-u n(r')r'dr' (20) 



The first current term H(^) arises in the normal Grad-Shafranov equation. The additional 
term in the GGS arises due to the electric field (see Eq. [17] ) acting on the moving struc- 
ture. This appears as an additional poloidal current, —uf n(r')r'dr', in the frame of the 
propagating solution, where the motionless ions are seen moving with velocity — ue z . This 
"current" can appear only where the electrons are moving. Furthermore, radial electric force 
or J x B force due to this "current" results in a pressure term — u f Q n(r')r' dr' x dH/d^/ 
in the GGS equation. It must be emphasized that in MHD the magnetic field is frozen into 
the plasma and cannot move relative to the ions, so there is no additional force. The GGS 
equations were first found in [l7] as exact solutions of the dynamic equations for the Hall 
magnetic field in an axisymmetric geometry, Eq. ( jTUl) and ( 1TTT) . 

IV. AXISYMMETRIC HALL MHD SOLUTIONS 

Here, we consider solutions of the GGS, that propagate down a cylindrical column. We 
have three sets of analytic solutions. Sec. ( 1IV Al) details the resistive nonlinear shock whose 



structure is described by a Burgers equation. Sec. fllVBI) describes a periodic solution to the 



GGS equation with linear current term H(^) oc Sec. (1IV CJ ) numerically estimates, with 
justification on the basis of the Chandrasekhar- Fermi theorem 29], nonlinear localized bubble 



solutions to the GGS integral. According to this theorem, a localized MHD equilibrium 



plasma magnetic structure could exist only in the presence of an external magnetic field or 
with fixed metal wall boundaries. 



A. Resistive Nonlinear Shock 



Again, assume a Hall plasma column, so that the magnetic fields are given by Eq. (TTT]) 
and ffHJ]) . with a const ant- density plasma within the column r < R and vacuum outside. 
Furthermore, let us assume only the toroidal field B exists with profile B(r, z, t) = B(z, t)r/R 
and B(—oo,t) = —B , where B > 0. Thus, we get the Burger's equation in normalized 
coordinates: 

dB _2BdB d 2 B 

~dt~lt~dz~ + m ~d^ (21) 
With this solution, in normalized coordinates: 

B*( r <R 



B(r,z,t) = { 2R \ V 1RV m / / (22) 
r > R 

In dimensional values the shock velocity is given by u = cB / (4nenR) and shock thickness 
5 = 2v ei /uj ce R. Note that the shock is moving in the positive z direction in our coordinate 
system. This asymmetry is typical for Hall phenomena. Furthermore, one may construct 



analytic solutions for any arbitrary initial configuration B(z,0) (see e.g. {4, |30| for the 
Hopf-Cole transformation that transforms Burger's equation into a heat diffusion equation). 
A plot of the toroidal magnetic structure of this resistive case is shown in figure [U with 
u = 100 and b = B /R = 100. 



B. Periodic GGS Solution 

Here, we consider solutions to the GGS equation (see Eq. [20] and [19]) with linear current 
terms, H(^f) = k^> . Note that this solution arises from a particular choice of current term, 
and does not arise from a linearization of the dynamic equations (Eq. [H] and [10]). Consider 
the density n = 1 a constant, and as before u > is also a constant. Then the GGS integrals 
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poioiaai current neia lines 
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FIG. 1: Poloidal current lines within the Burgers-like nonlinear solution. Here, electrons must 
return across the surface so, we see that a current sheet develops on the Hall plasma column, that 
distributes current back into the (attached) current generator. 



reduce to these relations inside the plasma column of radius R: 
dr \r dr J dz 2 2 



r— i -—I + 

B = kA— -ur 2 
2 



(23) 
(24) 



Outside the plasma column we have a vacuum, so that there are no currents and no magnetic 
field in the frame of electron MHD. Electron MHD implies that the magnetic field is carried 
only by the electrons. 

We consider a geometry in which there is a given radial variation in the three-component 
magnetic field, such that there is only a toroidal current sheet. This implies that B — > 



II 



at the boundary of the plasma column. We look for solutions that go as ~ cos/t (z — ut) in 
magnetic field. 



B 



uR 2 



ar.J\ (r^/k 2 — k 2 ) cos k [z — ut) + r < R 

r> R 

(20) 

«a Ji (rV k 2 — /c 2 ) cos k(z — ut) r < R 
r > i? 



(26) 



Where J n are Bessel functions, k 2 > k 2 , and the condition that Ji (-Ra/ k 2 — /c 2 ) = 0. The 
toroidal surface current on the plasma column is given by: 

1$ (R, z-ut) = «Vk 2 - k 2 ( J (r^k 2 - A; 2 ) - 

J 2 (rVk 2 - A; 2 ) ) cos k {z - ut) + — 

The solution shown in figures HI with these parameters: 

77 = 0.0 a =1.0 
k = 1.0 u — 100.0 

And k given by the first zero of J\[x): 



/3.83171V 



R ) 



C. Localized Hall Bubble Solutions 

Here, we consider solutions to the GGS equation that describe the evolution of possible 
2D nonlinear whistler-like waves and structures. Consider solutions with normalized density 
n = 1, propagating with velocity ue z , and the currents and toroidal magnetic fields in the 
GGS integral: 

H(V) = a^ 2 

1 1 ( 27 ) 

rB = H (\1>) ur 2 = a 1 ® 2 ur 2 

v ' 2 2 
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FIG. 2: Poloidal magnetic flux (and hence, poloidal magnetic field lines) for the GGS solution. 
Note that since the poloidal flux within the vacuum is constant, no poloidal magnetic fields exist 
within the vacuum. Note that the discontinuity in B z results in <j) current sheet along the Hall 
column surface. This magnetic configuration appears similar to the compact torus, where the 
structure is defined by the boundary conditions - i.e. metal walls. 
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In terms of the toroidal vector potential A and the toroidal magnetic field B the nonlinear 

soliton equation is given by, 

d 2 A IdA A 2 , n 2 43 n 
H - ur 2 A + 2r 2 aA i = 

(28) 



,2 1 

£> = «rA — -w 



2 

One can easily show that if we let A Q (r,z) be the solution to the differential equation, 
u — 1, a — 1: 

mr + —a" 2 - ~i ~ r A * + 2r A ° = 29 

Then the scaling for the toroidal vector potential A (r, u, a) in terms of A (r, z) and the 
form of the toroidal magnetic field from Eq. (|28~!) are given by: 

'a\i/2 



A (r, ^; /i, a) = (^) A (u l/4 r, u 1/4 z) (30) 

/q/\3/2 g I 

5 (r, z; a) — f — J rA (u lj r, u ' £j — -wr (31) 



Consider a volume to be a cylinder. The soliton structure consists of a three-component 
magnetic field (a poloidal field and a toroidal field) that is localized about the origin with 
a span much smaller than the dimensions of the cylindrical volume; thus, B r — > 0, B z — > 0, 
and B^ — > — ^ur at the boundaries of the cylinder. Employing the Chandrasekhar- Fermi 
theorem - using r • (J x B) force balance within the cylindrical volume - one can show, 

R rL 

/ rB 2 z dr dz = 

r=0 J z=—L /'QO\ 
R pL 

/ un{r)rH (\l/) dr dz 

r=0 Jz=-L 

Furthermore, since B z = —dA/dr, a localized solution (where the integral on the left hand 
side is finite) is possible provided that the right hand side > 0. To guarantee a localized 
solution we require that H(ty) > 0. The case analyzed in this section is H (\^) = \& 2 > 0. 

Note that the form of Eq. (128]) or ( 129]) is as if we consider a nonlinear 3D Schrodinger 
equation with the m = 1 azimuthal modes (i.e. solutions ~ e 1 ^). Here, to estimate the form 
of the soliton, we use the same variational method and trial functions, described below, for 
3D nonlinear optical solitons 



31 



321 ] . First, we derive the Lagrangian that describes the 
GGS equation with nonlinear current. Then we look for localized solutions with different 
test functions. One reasonable approximation for the test function in cylindrical coordi- 
nates is A(r, z) = U(r) sech \xz. A one-dimensional differential equation in r is derived by 

14 



averaging that Lagrangian with respect to z and searching for a localized one-node solution, 
presumed to be the main lowest-energy stable solution, for a given p. We solve for that p 
that minimizes the integral of the averaged Lagrangian over r. Another test function is in 
spherical coordinates, A(p,9) = U(p)sm9, where p = \/r 2 + z 2 and cos# = z/p. Again a 
one-dimensional equation in p, and its localized single-node lowest energy solution, is derived 
by averaging and minimizing the Lagrangian of the GGS. A more detailed description of the 
estimation method for both trial functions is provided in the Appendix. 

Plots of the dependence U(r), the cylindrical ansatz, and U(p), the spherical ansatz, 
are shown in Fig. [3J Contours detailing the shapes of the toroidal vector potential A(r, z) 
for both trial functions are shown in Fig. |H The cylindrical ansatz U(r) sech pz and the 
spherical ansatz U(p) sin 9 have approximate shapes in the radial coordinate. As in 311 ] . 
there is a substantial differences in the shape of the cylindrical and spherical trial functions 
of the soliton. The shape of the cylindrical and spherical trial functions look similar to their 
corresponding solutions of the optical solitons found in 31[. Both trial functions also satisfy 
the scaling A (r, z; p, a) = (a/u) 1 ^ 2 A (u 1//4 r, u l ^z) . 



V. DISCUSSION AND CONCLUSIONS 



Here we have explored possible structures that arise in a Hall plasma due to curvature 
rather than, say, density gradients. We have provided a general structure to continue further 
analysis of the dynamic Hall modes in Eq. (TTUT) and (TTTT) as well as the generalized Grad- 
Shafranov (GGS) integrals in Eq. ([19]) and fl20l) and propagating solutions with constant 
velocities in the z direction. 

Well-understood Hall plasma phenomena such as whistler drift modes were shown to 
exist in an homogenous cylindrical plasma column. Furthermore, an exact solution for the 
whistler drift mode, the Burgers nonlinear shock wave, was discussed in subsection II V A[ 
Periodic solutions were shown in subsection IIVBI from a simple analysis of the allowable 
solution for linear current in the GGS integrals. We find that the solution is constrained by 
the radial vacuum boundary of a Hall plasma column, because in Hall physics the magnetic 
field is carried solely by the electrons. Further work may involve a search for transitional 
solutions - those solutions intermediate between the periodic Grad-Shafranov solution and 
the collisional nonlinear shock. Possible solutions may be those with oscillatory substructure 
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FIG. 3: Plot of the radial distribution with the given ansatzes shown in the figure. For the 
cylindrical ansatzes, the z width parameter fj, is calculated to four decimal places. 

in the toroidal magnetic field at the shock with corresponding poloidal magnetic field in the 
case of small resistivity. 

We determined that there is a congruence between the GGS and the nonlinear Schrodinger 
equation with m — 1 azimuthal modes. We apply the variational method, which was devel- 
oped for optical fiber solitons, for our problem. We estimate soliton solutions in subsection 
IIV CI for a given nonlinear current H(^>) = a^ 2 and constant number density n. The 
Chandrasekhar-Fermi theorem implies that such a nonlinear current may allow for local- 
ized soliton solutions. Such a current distribution also allows for a simple scaling - that 
is, solutions with normalized propagation velocity m ^ 1 and current scalings a ^ 1, and 
u > l,cv > 1 have vector potential A(r, z;u,a) = (a/u) 1 Ao (u^r, v}l A z). Further re- 
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FIG. 4: Structure of the £7(r)sech^z ansatz (left) compared with the spherical U(p)sm9 ansatz 
(right). Although the finer features are obviously different, the width in the radial direction are 
roughly similar. 



search may apply this suggested method to other localized soliton-like solutions, as well as 
illuminate relationships between the magnetic fields and vertical velocities in propagating 
localized Hall structures. The Hall dynamic equations in Eq. ( jTUl) and ( ITT]) may be useful 
in analyzing the stability of the soliton solutions estimated here. 

We have also included some space plasmas in which the Hall regime may be important. 
In Earth's lower ionosphere, whistler drift modes due to density gradients (at the equator) 
or curvature (at the magnetic poles) may exist. In astrophysical plasmas, such as molecular 
clouds, protostellar disks, or protostellar jets, Hall phenomena may exist due to the presence 
of charged dust. 

Axisymmetric Hall phenomena may play an important role in the circumstellar disks 
and jets pictured by the Hubble Space Telescope around young stellar objects. The typical 
circumstellar disk has a radius of a few hundred AUs - a possible length scale in Hall phe- 
nomena. Hall shocks and other magnetic structures also transport of currents and magnetic 
fields through matter without compression or net mass motion. A comparison of the radio 
frequency Doppler studies of the bipolar jets around HH30 show that the bulk velocity of 



the matter is five times slower than the velocities of the observed jet substructures 



33). 
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Furthermore, radio polarization measurements of the jets emanating from the disks of GM 



Aurigae and DG Tauri imply largely toroidal magnetic fields [34]. which may be explainable 
as a nonlinear shock wave and structures rather than the traditional model of a supersonic 
fluid or MHD shock. 

Other explanations for the fast-moving, time-resolved substructures in the jets of T Tauri 
stars involve "bullets" of plasma that emanate from the disk 
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35 



8f. Furthermore, a variety 



36|]). Astrophysical MHD 



37 



381 ] ) may explain 



of MHD jet models are given in the literature (see, e.g. 
experiments of pressure-confined toroidal magnetic towers (see, e.g., 
features of these jets as specific MHD phenomena: higher intensity material within the jet 
as magnetic bubbles, and longer-timescale intermittency because kink and other azimuthal 
MHD instabilities collapse the jet structure. Here, due to the obvious presence of dust, 
we imply an alternative model, magnetic penetration into a Hall plasma, that may explain 
some of these observed structures. A companion astrophysical paper will expand upon the 
illuminating relationship between net Hall magnetic bubble velocity, radial dimension, and 
magnetic field strength to explain the migration of observed "bullets" seen in protostellar 
jets. 



Acknowledgments 



The author would like to thank L. Rudakov for his invaluable advice with regards to 
these new classes of curvature-based solution in Hall plasmas, and especially for introducing 
him to the generalized Grad-Shafranov equation. 



[1] A. S. Kingsep, Y. V. Mokhov, and K. V. Chukbar, Fiz. Plazmy 10, 854 (1984). 

[2] A. Fruchtman and L. Rudakov, Phys. Rev. Lett. 69, 2070 (1992). 

[3] J. D. Huba, I. M. Grossman, and P. F. Ottinger, Phys. Plasmas 1, 3444 (1994). 

[4] A. V. Gordeev, A. S. Kingsep, and L. I. Rudakov, Physics Reports 243, 214 (1994). 

[5] J. D. Huba, Phys. Plasmas 2, 2504 (1995). 

[6] M. A. Shay, J. F. Drake, R. E. Denton, and D. J. Biskamp, J. Geophys. Res. 103, 9165 (1998). 

[7] L. I. Rudakov and J. D. Huba, Phys. Rev. Lett. 89, 95002 (2002). 

[8] P. L. Pritchett, J. Geophys. Res. 106, 3783 (2001). 

18 



[9] J. D. Huba and L. I. Rudakov, Phys. Plasmas 9, 4435 (2002). 

[10] M. Wardle and C. Ng, Monthly Notices of the Royal Astronomical Society 303, 239 (1999). 

[11] S. Balbus and C. Terquem, Astrophys. J. 552, 235 (2001). 

[12] T. Sano and J. M. Stone, Astrophys. J. 570, 314 (2002). 

[13] T. Sano and J. M. Stone, Astrophys. J. 577, 534 (2002). 

[14] L. I. Rudakov, Phys. Scripta 2001, 158 (2001). 

[15] M. Salimullah and P. K. Shukla, Physics of Plasmas 6, 686 (1999). 

[16] M. Rosenberg and P. K. Shukla, Journal of Geophysical Research 105, 23135 (2000). 

[17] L. I. Rudakov, Phys. Scripta 2002, 58 (2002). 

[18] W. I. Newman, A. L. Newman, and R. V. E. Lovelace, Astrophys. J. 392, 622 (1992). 

[19] T. W. Harquist, W. Pilipp, and O. Havnes, Astrophysics and Space Science 246, 243 (1997). 

[20] B. T. Draine, W. G. Roberge, and A. Dalgarno, Astrophys. J. 264, 485 (1983). 

[21] S. Chapman, Nuovo Cimento 5, 1385 (1956). 

[22] S. Fromang, C. Terquem, and S. A. Balbus, Monthly Notices of the Royal Astronomical 

Society 329, 18 (2002). 

[23] P. C. Myers and A. A. Goodman, Astrophys. J. 326, L27 (1988). 

[24] M. Oppenheimer and A. Dalgarno, Astrophys. J. 192, 29 (1974). 

[25] L. Spitzer and M. G. Tomasko, Astrophys. J. 152, 971 (1968). 

[26] T. Umebayashi and T. Nakano, Monthly Notices of the Royal Astronomical Society 243, 236 
(1990). 

[27] E. G. Zweibel, Phys. Plasmas 6, 1725 (1999). 

[28] J. S. Mathis, W. Rumpl, and K. H. Nordsieck, Astrophys. J. 217, 425 (1977). 

[29] S. Chandrasekhar and E. Fermi, Astrophys. J. 118, 116C (1953). 

[30] S. I. Vainshtein, S. M. Chitre, and A. V. Olinto, Phys. Rev. E 61, 4422 (2000). 

[31] A. Desyatnikov, A. Maimistov, and B. Malomed, Phys. Rev. E 61, 3107 (2000). 

[32] A. Yakimenko, Y. Zaliznyak, and Y. Kivshar, Phys. Rev. E 71 (2005). 

[33] C. J. Lada and M. Fich, Astrophys. J. 459, 638 (1996). 

[34] M. Tamura, G. H. Hough, J. S. Greaves, J. -I. Morino, A. Chrisostomou, W. S. Holland, and 

M. Momose, Astrophys. J. 525, 832 (1999). 

[35] R. V. E. Lovelace, M. M. Romanova, and G. S. Bisnovatyi-Kogan, Monthly Notices of the 

Royal Astronomical Society 275, 244 (1992). 

19 



[36] G. V. Ustyugova, R. V. E. Lovelace, M. M. Romanova, H. Li, and S. A. Colgate, Astrophys. 
J. 541, L21 (2000). 

[37] A. Ciardi, S. V. Lebedev, C. Stehle, and T. Lery, in Journal de Physique IV (Proceedings) 
(LUTH, Observatoire de Paris, UMR 8102 CNRS, 92195 Meudon, France; The Blackett Lab- 
oratory, Imperial College, London SW7 2BW, UK, 2006), pp. 1043-1045. 

[38] A. Ciardi, S. V. Lebedev, A. Frank, F. Suzuki- Vidal, G. N. Hall, S. N. Bland, A. Harvey- 
Thompson, E. G. Blackman, and M. Camenzind, Astrophys. J. 691, L147 (2009). 



Appendix: Variational Method For Trial Functions of GGS 

Here, we demonstrate localized 3D soliton solutions of the GGS equation with nonlinear 
current term, with the cylindrical coordinate and spherical coordinate trial functions. This 
variational method and trial functions were borrowed from jsi|. Furthermore, the integration 
of the one-dimensional functions was performed using a fourth-order Runge-Kutta method 
from r = to r = 10 for the cylindrical ansatz and p = to p = 10 for the spherical ansatz. 



1. Cylindrical Coordinate Ansatz 

The Lagrangian of the differential expression (1291) is given by, 

L (d r A , d z A , A o] r, z) = K (d r A ) 2 + l -r (d z A ) 2 + 4f + ±r 3 A 2 - ±r 3 A 4 (A.1) 
With this ansatz, 

A (r, z) =U (r) sech (pz) , (A. 2) 
we get the averaged Lagrangian C±, where U' = dU/dr: 

d = T L (d r A , d z A , A ; r, z) dz = - (f-) * + ^U 2 + — + -U 2 - ^ 4 3) 
J-oo p \dr J 3 pr p 6p 

If we vary C with respect to U and dU jdr we then have a differential equation in r: 

d 2 U ldU U 1 2tt 2tt 4 2rr3 n 

— + — -p 2 U - r 2 U + -r 2 U 3 = (A.4) 

dr z r dr r z 3 3 

The solutions that we look for are given by the 
limU(r) = U'(0)r 

r— >0 

/ i \ (A-5) 
lim U(r) = r 1 exp I — -r 2 

r— >oo \ 2 



20 



We employ the shooting method, with free parameter U'(0), by integrating out Eq. (1A.4j) 
from r = to some sufficienly large r = ru m 3> w, where w is the width of the structure. 
Furthermore, since we postulate the lowest-energy solutions are probably dominant, we look 
for solutions that have a single node in U(r). 

For each /i there is a unique single-node solution U(r;fi) satisfying Eq. (lA.4j) . To solve 
for this problem, we require that the averaged Lagrangian Eq. (1A.3j) . averaged over r, must 
be extremized with respect to \i: 

r°° dd 

Jo dfx 



dr = 



(A.6) 



If we use these functions: 



ei 0) 
es {fJ-) 



rU (r; /i) 2 dr 

i 

r 3 U (r; /i) 2 dr 

i 

r 3 U (r; /i) 4 dr 



(A.7) 



We then get the relation 

duV r 2 

'o \ \ dr J r 
However, using flA.4j) one can derive this relation, 



— ) zj ) dr = ^/i 2 ei (//) - c 2 (/') di). 



2 
3 ( 



(A.8) 



— 

dr J r 2 



A. ( rU ^L) _ I^ r f/ 2 - r 3 U 2 + -r 3 U 4 . 
dr \ dr ) 3 3 



(A.9) 



Substituting Eq. (|A.9[) into the integral on the left hand side of Eq. flA.8p . and noting that 
the first right hand side term in Eq. flA.9j) integrates to zero because it is a surface term, we 
then get the residual function for /i : 



(A. 10) 



Our estimates of /x (to four decimal places) and Uq that satify flA.lOp are given below: 
u w 2.1852 

(A.ll) 

Z7'(0) w 1.0115 
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2. Spherical Coordinate Ansatz 



For the spherical coordinate ansatz, we transform (|29|) into spherical coordinates, with 
p = y/r 2 + z 2 and cos 9 — zj y/r 2 + z 2 : 



1_8_ ( 2 cMo\ 1 8 ( . 8Aq 
p 2 dp V dp J p 2 sin 9 89 V m 89 

4 " p 2 sin 2 9A + 2p 2 sin 2 9A 3 = 



(A.12) 



p 2 sin 2 9 
Which has the Lagrangian, 



2 r V dp 



2 



L(8 p A ,8 e A ,A ;p,9) = \p 2 sm9 ) i 



± sin ( ^ 2 + _4L + I p 4 sin s M 2 _ 1 4 gin 3 M 4 
2 I 90 / 2sin# 2^ 2^ 



(A.13) 



With the ansatz: 

A o (p,0) = U(p)sin9 (AAA) 
We get the averaged Lagrangian: 

C 2 = £l (8 p A , 8 e A , A 0] p, 9) d9 = \p 2 {^\ * + -U 2 + ^p'U 2 - f/^-^) 

That results in this differential expression: 

d 2 U + 2dU _2U_ 4 2(J + 48^ 2[/3 _ Q 
dp 2 p dp p 2 5 35 

We shoot with U(0) = and free parameter U'(0) such that for sufficiently large p, we get 
U(p) —> ap^ 1 exp (— p 2 /5), where a is some constant. U'(0) is given by, 

U'(0) « f.5897 (A.17) 
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